Non-Fermi Liquid Behavior and Double-Exchange Physics 
in Orbital-Selective Mott Systems 
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We study a multi-band Hubbard model in its orbital selective Mott phase, in which localized elec- 
trons in a narrow band coexist with itinerant electrons in a wide band. The low-energy physics of 
this phase is shown to be given by a generalized double-exchange model. The high-temperature dis- 
ordered phase thus differs from a Fermi liquid, and displays a finite scattering rate of the conduction 
electrons at the Fermi level, which depends continuously on the spin anisotropy. 

PACS numbers: 71.27.+a, 71.30.+h, 71.10.-w 



The Mott phenomenon, localization of electrons by 
their mutual interactions, is a key phenomenon in corre- 
lated electron materials. Recently, there has been great 
interest in the possibility of an "orbital-selective Mott 
phase" (OSMP), in which a Mott transition occurs within 
a restricted subset of quasi-degenerate orbitals, so that a 
narrow band of localised electrons coexists with a wider 
band of itinerant electrons. Originally proposed for 
the ruthenate compound Ca^S^-zRuO^, this possibility 
has been intensively discussed since then |2|, y, y, |3, |fj|, |7| • 
In this letter, we clarify the nature of such an orbital- 
selective Mott phase. We show that for finite Hund's 
coupling the orbital-selective Mott transition is accom- 
panied by a breakdown of Fermi liquid theory, which can 
be understood from a mapping onto an effective double- 
exchange model at low energy. In the high-temperature 
disordered phase, the electrons in the conduction band 
acquire a finite lifetime, due to the scattering on the 
localised spins. This behavior is eventually cut-off by 
long-range ordering at low temperature. Detailed calcu- 
lations within dynamical mean-field theory (DMFT) are 
performed in order to demonstrate these properties. 

We consider the two-orbital Hubbard model for a 
broad and a narrow band, described by creation oper- 
respectively, with hamiltonian: 
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In this expression, U and U' are Coulomb interaction pa- 
rameters and J is the Hund's rule coupling. We have 
used the notation — c\ a Ci,y — 1/2 (similarly, nf a ) 



and written the hamiltonian in the manifestly particle- 
hole symmetric case corresponding to two electrons per 
site (n c = n d = 1), on which we shall focus in the 
following. The spin- flip and pair-hopping terms of the 
hamiltonian read: Hfi ip = 



"^2i[ c iT Ci l4i c ^T + n - c ] and 
■f h.c] . A variable parameter 



a has been introduced in front of these terms, for fur- 
ther use. For a cubic environment and unbroken spin 
symmetry, the relations U' = U — 2 J and a = 1 hold. 
This model has been the subject of several recent stud- 
ies using dynamical mean-field theory (or equivalently, in 
the limit of infinite lattice connectivity) . For anisotropic 
interactions (U' < U and finite J) an orbital selective 
Mott phase is found even at moderate values of the ratio 
of the two bandwidths Q . As shown recently 0, , for 
large differences in the bandwidths the OSMP can be in- 
duced even in the case of SU (4)-symmetric interactions 
(U = U',J = 0). 

First, we consider this model for a = 0, i.e with 
density-density interactions only, and make a compar- 
ative study of the nature of the paramagnetic OSMP 
in the absence (J = 0) or in the presence (J ^ 0) of 
the Hund's coupling. The model is solved using DMFT, 
for two semi-elliptical densities of states (infinite connec- 
tivity Bethe lattice) p 0c (e) = 2[1 - (e/Dj 2 ] 1 / 2 /^) 
(similarly pod), with a bandwidth ratio Dd/D c = 0.1. 
We use D c = 1 as our unit of energy, choose U = 0.8, 
U' = U — 2J, and compare J = to J = 0.2. In both 
cases, the self-energy of the narrow band (not shown) 
diverges at low frequency, while the self-energy of the 
broad band does not, indicating that the narrow band 
is Mott-localized while the broad band is not (OSMP 
regime). However, the low- frequency behavior of the 
self-energy is very different in each shown on 

Fig. n which displays our DMFT results obtained with 
the QMC algorithm of Hirsch and Fye|g. It is seen that, 
for J = 0, the self-energy extrapolates to zero at low fre- 
quency (ImS c (iu re ) ~ (1— 1/Z c )u>„+- ■ ■) in a Fermi-liquid 
manner, while for J =/= 0, it extrapolates to a finite value 
ImE c (zO + ) = T c 7^ 0. Hence, a finite lifetime is found at 
the Fermi level for J ^ 0, so that no well-defined Fermi- 
liquid quasiparticles exist. Correspondingly, this implies 
a violation of the Luttinger sum-rule p c (0) — Pco(0) for 
the spectral function p c {u>) = —hsxG c {u> + iO + )/ir. The 
local Green's function being related to the self-energy by: 
G c {iu> n ) = J dep c o(e)/(iu/ n —H c (iu) n ) — e), anon-zero life- 
time implies: p c (0) = J dep c0 (e)T c /[Tr(e 2 +T 2 C )} < p c0 (Q). 
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On Fig. [21 we display the spectral functions of the nar- 
row and broad bands, obtained by a maximum-entropy 
continuation of the QMC data. Pinning at the Luttinger 
theorem value p c (0) = p c o(0) = 2/(nD c ) holds for J = 
but not for J / 0. This is consistent with Fig. 3 in 0- 
Note that the spectral function of the (localised) narrow 
band has a sharp Mott gap in the latter case, but dis- 
plays rather a pseudogap in the former (as previously 
discussed in 0,0]). The Green's functions obtained for 
J = 0.2 using QMC and using an exact diagonalisation 
(ED) solution of the DMFT equations (inset of Fig. [TJ 
are in perfect agreement. 
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FIG. 1: Broad band self-energy ImE c (i(-On) for different tem- 
peratures ((3D C = 40, 60, 80, 120 correspond to circles, empty 
and filled squares and diamonds respectively). In the J = 
case (online: black) Im£ c (iaj) extrapolates linearly to zero 
(Fermi liquid behavior), whereas for J = 0.2 (online: red) a 
finite lifetime is found. Inset: ImGd (iui„ ) and ImG c (iui„) for 
U=0.8 D c and J=0.2 D c . QMC results (symbols) correspond 
to finite temperatures ((3D C — 40, 60, 80, 120 given by circles, 
squares, triangles up and down respectively), while the ED 
data (lines) are at zero temperature. 

In order to understand the origin of this finite-lifetime, 
iron- Fermi liquid behavior, we now derive an effective 
hamiltonian valid at low-energy in the orbital-selective 
regime. There, it is expected to be legitimate to ne- 
glect entirely the hopping term of the narrow band (which 
renormalizes to zero at low-energy). Furthermore, since 
we are interested in low-energy physics, the states with 
a hole or a double-occupancy in the narrow orbital can 
be eliminated. Setting t d — and n d = ±1 in @, a 
straightforward calculation leads to the effective hamil- 
tonian: 

#cff = tC ijA„ c 3« + U H "it"! - 

ija i 

-2JJ2 [SUSto + « + S v ld St)] (2) 

i 

In this expression, S^ ,z are the components of the 
spin- 1/2 operator acting on the 2 local moment states 



FIG. 2: Spectral functions for J = and J = 0.2 (U = 
0.8D c ,f3D c = 120) calculated from QMC. The narrow band 
Pd(u) (dashed line) displays a gap or pseudogap in both cases. 
For J — 0, p c {uj) (solid line) coincides at ui = with the non- 
interacting density of states (dashed-dotted line), as a conse- 
quence of Luttinger theorem, while this condition is clearly 
violated at J / (p c (0) < poc(0)). 



I T)d, I i)d at each site, while Sf c = l/2(c] T c iT - c^c^), 
Sf c + iS^ c = S£. = cLcjj, are the components of the itin- 
erant electron spin operator. Note that the inter-orbital 
Coulomb interaction U' entirely disappears from this ef- 
fective hamiltonian (as well as the pair-hopping term, 
which cannot act in the low-energy subspace of sites with 
single-occupancy of the d-band). Hence the low-energy 
effective hamiltonian in the orbital-selective phase is a 
ferromagnetic Kondo lattice, with an additional Hubbard 
interaction in the (broad) itinerant band. The orbital- 
selective Mott localisation of a single orbital component 
generates an effective double-exchange model in which 
the itinerant orbital interacts with the "core-spin" of the 
localised orbital through the Hund's coupling. This ob- 
servation has two important consequences, to be explored 
in more detail below: (i) it explains the observed non- 
Fermi liquid behavior, and its dependence on the cou- 
pling to spin-flip terms a and (ii) it suggests that ulti- 
mately the orbital-selective phase undergoes long-range 
ordering as temperature is lowered. 

When the interactions are of the density-density form 
(a = 0), the itinerant orbital couples to the core-spin 
through Ising terms only. The effective hamiltonian @ 
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then essentially reduces to a Falicov-Kimball model for 
each spin species, in which the conduction electrons in- 
teract with a given configuration of the Ising variables 
Sf z , which are conserved quantities. In addition, the two 
spin components of the itinerant band are coupled by 
the intra-orbital Hubbard interaction. In order to solve 
this model, one must consider a given configuration of 
the Ising spins, solve the one-band model of conduction 
electrons in that configuration, and average over all Ising 
spin configurations with appropriate Boltzmann weight. 
This process can be achieved within DMFT (infinite co- 
ordination) by mapping the problem onto an effective 
single-site action: 

Sg° = - jdr J dr'^2(ci(r)g^(r-r')c a (r') 

± / drpn^ni- JS z {n^ -nj] (3) 



subject to the self-consistency condition Qq (iu> n ) — 
iuj n — tlG c {iuj n ). Because S z is a conserved quantity, 
each sector S z — ±1/2 can be considered independently. 
In each sector, the effective action @ is that of an ef- 
fective Anderson impurity model in a local magnetic 
field ±J. The itinerant electron Green's function is ob- 
tained as: G c — l/2[Gj+ + Gf_], where G a ± denotes 
the Green's function — (Tc a c a )± in the effective action 
[ffil with S z = ±1/2, and the obvious symmetry prop- 
erty G a + = G S - holds. In writing these expressions, a 
paramagnetic state has been assumed and particle-hole 
symmetry corresponding to half-filling has been used. In 
Fig-03 we compare the Green's function obtained by solv- 
ing this effective model to the full solution of the original 
model in the OSMP. The two results are indistinguishable 
for all practical purposes, hence firmly establishing the 
validity of our low-energy effective model in this phase. 
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FIG. 3: Green's function in imaginary time for U=0.8D C , 
J=0.2D C , j3 — 80, computed from QMC. The red curves give 
the results for the two-band model; the effective one-band 
model Green's function (empty black circles) fit the wide band 
of the two-band model perfectly. The green curves represent 
the auxiliary Green's functions of the S z = ±1/2 sectors of 
the effective 1-band model. 



For a given orientation of the core spin (say, S z — 
±1/2), the Green's functions of the itinerant electrons 
for each spin component correspond to those of an impu- 
rity model in a local magnetic field jJJJ, an d hence read: 
G T+ (iw„) _1 = G^ 1 ± J/2 - S T+ (icj„), Gx + (icj„) _1 = 
Qq 1 — J/2 — E|_|_(icj ra ), with the symmetry relation (at 
half- filling): Sf+(icJ n ) = — iui n ). In these expres- 

sions, E CT _|_ is the self-energy originating from the intra- 
orbital Hubbard interaction. In the itinerant (OSMP) 
phase for the conduction electrons, these spin-resolved 
self-energies for a given core spin orientation do have a 
Fermi-liquid form at low energy: E| + (iw n ) = En ± (1 — 

l/Z)ioj n H , with En a real quantity [ijj- J = J — 2En 

can be viewed as the effective exchange coupling renor- 
malised by the on-site Hubbard interaction. Inserting 
these expressions into G c = l/2[Gj+ ± Gj,+], and tak- 
ing into account the self-consistency condition Qq 1 = 
iuj n — t 2 c G c allows us to obtain the low-frequency form 
of the full self-energy E c in the form: 

E c (w ± i0+) = -iT c ± (1 - l/Z eS )oj + ■■■ 

r — J 7.- 



2(4t 2 -J 2 )V 2 



z -l (8^-^)(2t^J- i ) 

(4*2-J2)2 1 ' 



Hence, a finite lifetime at zero energy emerges, due to 
the scattering of the itinerant electrons onto the core 
spin associated with the localised orbital. Because the 
core S z is a conserved quantity in the absence of spin- 
flip (a — 0), the physics is closely related to that of the 
Falicov-Kimball model, in which this kind of non-Fermi 
liquid behavior was first emphasized by Si et al. , and 
later by Furukawa |12J for double-exchange models in the 
context of manganites (see also |l3[ E3 . Ila llrJ | ). 

An obvious question is whether this finite lifetime is a 
robust feature when spin-flip terms (a ^ 0) are reintro- 
duced. In that case, we have to deal with a ferromag- 
netic Kondo lattice model, with full 5/7(2) symmetry at 
a = 1 and Ising anisotropy at a < 1. Within DMFT, this 
model maps onto a ferromagnetic Kondo impurity prob- 
lem, with a self-consistent conduction electron bath. The 
renormalisation flow of the ferromagnetic Kondo problem 
[l7l | is such that the spin-flip component of the interac- 
tion is irrelevant as long as a < 1, with a continuous 
family of infra-red fixed points corresponding to X-ray 
edge physics. Hence, we expect a finite lifetime to sur- 
vive as long as the SU(2) symmetry is broken (i.e for 
all a < 1). We have solved the DMFT equations for 
the original 2-orbital model in the presence of spin-flip 
(and pair-hopping) terms for a range of values of a, in 
the orbital-selective phase. The results are displayed in 
FigEl an d indeed fully support this expectation, with a 
violation of Luttinger's sum rule gradually decreasing as 
a is increased. In the fully symmetric case a = 1, the 
lifetime does extrapolate to zero (FigBJ as w -* 0, but 
conventional Fermi liquid behavior is not recovered. In- 
deed, single-impurity Kondo scaling suggests a singular 
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behavior Im£ c ~ — l/(lnw) 2 , consistent with the fit of 
our numerical results displayed in Fig. [3] 
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FIG. 4: Low-energy part of ImE c (iw), obtained by ED at 
T = for U = 1.2D C , J = 0.15(7, and different values of a 
(0, 0.3, 0.5, 0.7, 1 from bottom to top). It extrapolates to zero 
at u) — only in the case of fully symmetric interactions (a = 
1, red solid curve). Inset: Zero-frequency extrapolation of the 
Green's function as a function of J and for different values of 
a — 1, 0.7, 0.5, (from bottom to top). In the case of a fully 
SU(2) symmetric interaction vertex (a = 1) the Luttinger 
pinning condition ImG(0) = —2/D c is fulfilled. As soon as 
this symmetry is broken (a < 1), the pinning condition is 
lost. 

The calculations described above are restricted to 
the paramagnetic phase, spin symmetry being enforced 
on the Green's functions. Clearly, at low tempera- 
tures, a phase with long-range order is more favorable. 
The competition between ferromagnetism, antiferromag- 
netism and phase separation has been studied in detail 
for double exch ang e models relevant to manganites (see 
e -g 113 Q Fewer studies have taken into ac- 

count the effect of an on-site Hubbard interaction, which 



increases the tendency to ferromagnetism [l9j |. Such a 
detailed study for the present model is beyond the scope 
of this paper. We have however performed some calcu- 
lations without spin- symmetrisation and do find a fer- 
romagnetic solution to be stable below a temperature of 
order D c /50. The inset of Fig.[S]demonstrates how a fer- 
romagnetic state, by suppressing spin disorder, restores a 
scattering rate which decreases as temperature is lowered 
and as the magnetic polarisation increases. 

In conclusion, we have shown that the low-energy 
physics of the orbital-selective Mott phase is given by a 
double exchange model (with an additional Hubbard re- 
pulsion). As a result, in the high-temperature disordered 
phase a non-Fermi liquid state with a finite lifetime of 
the conduction electrons is found. At lower temperature, 
long-range spin ordering and possibly phase separation 
sets in. By introducing a finite crystal field splitting, 
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FIG. 5: Low-energy behavior of ImE c (ioj) obtained by ED for 
U = 1.2D C and J = 0.1517 (squares), J = 0.2517 (circles) in 
the symmetric case a=l. Solid lines are low- frequency fits of 
the data to the function F(u) = — a/[(\ncu) 2 +b\nu>+c]. Inset: 
ImE c (w) for U = 0.8D C , J — 0.2D C , a — 0, computed from 
QMC without spin symmetrisation. Empty and solid circles, 
squares and triangles (up and down) correspond to j3D c = 
10, 20, 30, 40, 50, 120. The last two temperatures correspond 
to a ferromagnetic state. 



it is clearly seen that the double- exchange physics and 
that of the OSMP are continuously connected, the for- 
mer corresponding to a large crystal field splitting and 
a core spin being formed in the lower crystal field level, 
while the latter corresponds to a small crystal field and 
Mott localisation of the narrow band. Conversely, the 
OSMP is a useful way of simulating the double exchange 
model with fully quantum spins using DMFT techniques. 
Finally, it is tempting to speculate that a selective locali- 
sation of electrons in momentum space may lead to short 
lifetimes over part of the Fermi surface, as observed in 
the cuprates. 

Calculations were performed at IDPJS Orsay (project 
number 051393). We thank A. Millis for discussions. 
While this work was being written, independent work 
by A. Liebsch suggested the OSMP to be a bad metal 
phase with non- Fermi liquid behavior |2fjj. 
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